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If pj = f>j + Apj (j=l, M) y and y/ = y/ 0 + Ay/ , then equation 2 can be written, 
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Abstract: 

This work describes a new approach to multi-point "Dixon" fat-water separation that is 
amenable to pulse sequences that require short echo time increments, such as steady state 
free precession (SSFP) and fast spin echo (FSE) imaging. Using an iterative linear least 
squares method that decomposes water and fat images from source images acquired at short 
echo time increments, images with high SNR and uniform separation of water and fat are 
obtained. This algorithm extends to multi-coil reconstruction with minimal additional 
complexity. Examples of single coil and multi-coil fat-water decompositions are shown 
from source images acquired at both 1.5T and 3.0T. Examples in the knee, ankle, pelvis, 
abdomen and heart are shown, using FSE, SSFP and spoiled gradient echo pulse sequences. 
The algorithm was applied to systems with multiple chemical species, and an example of 
water-fat-silicone separation is shown. An analysis of the noise performance of this method 
is described and methods for improving noise performance through multi-coil acquisition 
and field map smoothing are discussed. 



Key Words: fat suppression, musculoskeletal imaging, cardiac imaging, magnetic 
resonance imaging, steady-state free precession, fast spin echo, phased array coils, silicone 
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Introduction: 

Steady-state free precession (SSFP) is a rapid, short TR imaging technique with 
specific advantages compared to short TR gradient echo techniques, including a high 
signal-to-noise ratio (SNR) and favorable contrast behavior, especially for visualizing fluid, 
because its contrast depends upon both Ti and T 2 [1-3]. However, utilization of SSFP has 
been limited by the fact that fluid and fat both appear bright on SSFP images. This 
characteristic of SSFP may cause abnormalities to appear similar to normal fat and, 
thereby, obscure underlying pathology. 

The application of "Dixon" fat-water separation method to SSFP imaging is an 
approach that could potentially provide homogeneous and reliable separation of fat and 
water from SSFP images [4, 5]. Current methods for SSFP fat suppression include 
fluctuating equilibrium magnetic resonance (FEMR), linear combination SSFP, and fat 
suppressed SSFP; however, all are sensitive to field heterogeneities [6-8]. The combination 
of "Dixon" methods with SSFP has been challenging for several reasons. First, SSFP 
requires short repetition times (TR) to prevent image degradation from field heterogeneities 
[1, 2], a constraint that limits TE increments to values that are smaller than those 
traditionally used in three point "Dixon" methods [5]. In addition, resonant frequency 
offsets from chemical shift and field heterogeneities produce additional phase shifts unique 
to SSFP that are problematic for "Dixon" fat-water decomposition techniques [2, 6]. 

Application of "Dixon" imaging to fast spin-echo (FSE) sequences has also been 
limited because the acquisition of echoes at different time shifts with respect to the spin- 
echo increases the spacing between successive refocusing pulses (echo spacing) [9]. 
Increasing the echo spacing reduces the number of echoes that can be collected in a time 
that maintains acceptable blurring from T 2 decay [10], offsetting the scan time benefits of 
FSE. A fat-water separation method that permitted shorter time increments would reduce 
the time between refocusing pulses and be beneficial to fast spin-echo imaging. 

Dixon fat-water decomposition techniques have historically been limited to single 
coil acquisitions because decomposition algorithms require phase unwrapping algorithms 
that are problematic when multiple surface coils are involved. Fat-water decomposition 
with multi-coil reception in combination with phase unwrapping algorithms has been 
described recently [11]. 
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In this work, we describe a new method for implementing "Dixon" fat-water separation 
utilizing an iterative least-squares method that reconstructs data acquired at short echo time 
increments, yielding images with high SNR and uniform separation of fat signal from water 
signal. This algorithm extends naturally to multi-coil reconstruction with minimal 
additional complexity and single and multi-coil decompositions derived from images 
obtained at both 1.5T and 3.0T are shown. Examples of fat-water separation in the knee, 
ankle, pelvis and heart are shown with SSFP, FSE and other pulse sequences. The 
algorithm is also applied to separation of multiple chemical species and examples of water- 
fat-silicone separation are shown. An analysis of the noise performance of this method is 
provided and methods for improving noise performance through field map smoothing are 
discussed as well as the effects of errors in presumed chemical shift are discussed. 

Theory: 

The application of traditional three-point "Dixon" fat-water separation to SSFP and 
FSE has been challenging. Three point methods described previously by Glover [5] 
describe the special case of phase shifts of 0, and In which correspond to echo time 
increments of 0, 2.2ms and 4.4ms at 1.5T and 0, 1.1ms and 2.2ms at 3T. Unfortunately, 
such echo time increments cause significant lengthening of the minimum TR. When using 
SSFP, increases in TR can lead to severe image degradation from banding artifacts caused 
by off-resonance field heterogeneities and chemical shift. When applied to FSE, this 
method increases echo spacing by 4.4ms at 1.5T, limiting the maximum echo train length 
in order to prevent blurring from T 2 decay, and also the number of slices per unit scan time. 
A recently described modified three point method that acquires images at phase shifts of 0, 
n/2, and /r, acquiring images at 0, 1.1, 2.2ms at 1.5T has been applied to FSE and partially 
alleviates this problem and is restricted to a specific echo time increment [12]. 

An iterative linear least-squares approach was formulated and a generalized 
algorithm with arbitrary echo times and multiple chemical species was developed and is 
described below. This is followed by a description of its extension to multi-coil 
applications and an algorithm that summarizes our approach to decomposition of each 
chemical species. Finally, a noise analysis of this algorithm is provided. 
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Signal Model 

Consider the signal in an image from a pixel containing M species each with 
chemical shift (Hz), Afj (j=l, M) located at position, r, acquired at an echo time, U 



M 



where pj is the intensity of the f species, and is in general, a complex term with its own 
magnitude, \p, I and phase, fa and y/ is the local magnetic field heterogeneity (Hz). If 
measurements are made at discrete echo times, t„ (n=l, .... N), then, 

representing the signal in a pixel located at position r, at echo time, t„. For FSE, t n is the 
shift from the center of the spin echo that occurs at <„=0. Eq. 2 contains M complex 
unknowns (pj ,j=l,...M)) and one scalar unknown O), for a total of 2M+1 unknowns. 
Each image contributes a real and imaginary measurement, constituting two measurements 
per time point, t n . Therefore, in general, M+l or more images are required to determine the 
system and separate all chemical species. For example, in a system with fat and water, 
M=2, and at least three or more images are required to separate fat and water. 

Least Squares Estimation of Fat and Water Images 

If an initial estimate of the field map, i// Q is known, then eq. 2 can be rewritten, 

M 

Eq. 3 is a linear system of complex equations that can be split into real (5* ) and imaginary 
(s'„) parts, 

where p* and p) are the real and imaginary components of the / species, 
c Jn =cos(2«A//J and d Jn = sin(2«A//„). Eq. 4 forms a set of linear equations that is 
amenable to linear least squares fitting to decompose estimates of each chemical species. A 
detailed description is provided in appendix A. 
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The initial estimate of the field map (y/ 0 ) is then refined further by defining error 
terms: y/=y/ 0 + &y/ 9 p*=p*+&Pj, and p[ = p] + Hp) Q'=1....M). Inserting these 
expressions into eq. 2, estimates of Ay/ , Ap*andApj can be calculated in the least- 
squares sense. Details are described in appendix B. Mathematics particular to the special 
case for systems that contain only water and fat are described in appendix C. 

Algorithm for Single Coil Acquisitions 

Using the above equations and those in the appendices, the following algorithm 
summarizes the method used to determine the least-squares estimates of water images and 

fat images for each pixel: 

o Estimate the signal from each chemical species using eq. A.3 and initial 

guess for field map, y/ 0 . A useful initial guess for y/ Q is zero (Hz) 
o Calculate error to field map, Ay/, using eq. B.7 
o Recalculate i// = y/ 0 + Ay/ 

o Recalculate s n (eqs. 3 and 4) with the new estimate of y/ 

o Repeat preceding three steps until Ay/ is small (eg. <1 Hz) 

o Spatially filter (smooth) final field map, y/, with low pass filter 

o Recalculate final estimate of each chemical species images with eq. A.3 

Filtering of the final field map is used to improve noise performance and is discussed 
below. 

Multi-Coil Acquisition and Reconstruction Algorithm 

A multi-coil acquisition with P elements collects P independent images, all with a 
different relative phase offset. By using the algorithm described above, P images of each 
chemical species are generated, as well as P field heterogeneity maps, which should not 
depend on coil-dependent phase shifts. The field maps are then combined by weighting the 
contribution from each coil by the square of the magnitude of the image contributed by that 
coil, in a manner similar to that performed for standard multi-coil image combination [13- 
15]. Specifically, for each pixel the combined field map is calculated as, 



6 



Vc = 




(5) 



Using the combined field map s n is recalculated with eq. 4 and new estimates of the P 
images of the different chemical species are calculated from eq. A.3. Finally, the P images 
for each chemical species are combined using a commonly used multi-coil reconstruction, 
previously described [13-15]. 

A summary of the final chemical species decomposition is provided below, 
accounting for multi-coil acquisitions, and smoothing of the final field map, 
• For each coil 

o For each pixel 

■ Estimate each chemical species assuming initial guess for field map, 
y/ Q , from eq. A.3 



• If this is a multi-coil acquisition, combine P field maps into combined field map 
using eq. 5 

• Spatially filter (smooth) field map with low-pass filter 

• Recalculate final images for each chemical species with eq. A.3 for each coil 

• For multi-coil acquisitions, combine final images for each chemical species from 
each coil using standard multi-coil combination techniques. 

Noise Considerations 

In the presence of noise with variance o*, the covariance matrix of the estimate of 

p fromeq. A.3 is [16], 



if it assumed that the field heterogeneity map, y/, is known. The diagonal elements of 
cov(p) represents the Cramer-Rao bound, or lowest possible limit of the error variance of 



Calculate error to field map, Ay/, from eq. B.7 



Recalculate y/ = y/ 0 + Ay/ 

Recalculate S n (eqs. 3 and 4) with the new estimate of if/ 
Repeat process until Ay/ is small. 



cov(p) = cr 2 (A T A)- 1 



(6) 
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the ht h species. In the specific case of water and fat, with equal echo time spacing of 
increment At, such that t n =t 0 +nAt, with «=0,1, ...N-l, the variance of fat and water 
estimates are equal and are determined from the diagonal terms of the covariance matrix 
calculated from eq. 6, 

Var( ^ /2 (7) 
sin 2 9 

where 9 = n Af^At , and it has been assumed that real and imaginary parts of the signal are 
uncorrected and have equal variance. As described by Glover [5], the effective number of 
signal averages (NSA) is, 

var(p) 

and is a helpful way to describe the effect of echo time on noise behavior of a multi-point 
fat-water separation. Therefore, for evenly spaced echo times, the effective number of 
signal averages is, 

1 sin 2 NO ( q\ 

NSA = N r— W 

N sin 2 9 

representing the upper limit of effective signal averaging and noise performance for any Ap- 
point fat-water estimation technique. The optimal echo time spacing that maximizes NSA 
occurs when N9 = 7t, such that At = fatf - which is an intuitive result reflecting the 

fact that optimal sampling should occur when the phase differences between fat and water 
are evenly distributed around the unit circle. At this echo spacing, NSA=N, demonstrating 
that for images acquired at the optimal echo spacing, the SNR of the calculated water and 
fat images is equivalent to average of N source images. Therefore, multi-point "Dixon" 
decomposition performed at the optimal echo time spacing is efficient from a signal to 

noise ratio perspective. 

Figure la plots eq. 9 for N=3 and N=4 at 1.5T, representing the upper limits of 
noise performance for a multi-point fat water decomposition, assuming that the relative 
chemical shift of fat and water is -3.5ppm (A/^ =-220 Hz at 1.5T). For example, the 
optimal echo spacing at 1.5T for a three equally spaced echoes is 1.5ms, and 1.1ms for a 
four echo acquisition. 
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In the case of three species, the calculated values of NSA are not equal for each of 
the three species. For such a system, the NSA for each of the three species is: 



NSA. 



_ ::i kABC r N{4 2 2 + B*)) (10a) 



N -C 



NSA,= N M BC ^* C '1> (10b) 
5 N -A 

sin(^VA0 = sinW^AQ = sin(M/ 23 ^A0 and ^ 
where ^- • sin^A/nAO sin(^A/ 23 Af) 

A/ 23 are relative chemical shifts of species 1 to 2, species 1 to 3, and species 2 to 3, 
respectively. Fig. lb plots eqs. lOa-c assuming the chemical shifts of fat and silicone 
relative to water are -3.5ppm (-220Hz at 1.5T) and -4.9ppm (-310Hz at 1.5T), respectively 
[17]. From this figure, the increment in echo time that optimizes the noise performance for 
silicone, fat and water is approximately 2.4ms at 1.5T. 

Materials and Methods: 

Images of human volunteers were obtained on a 1.5T GE CV/i scanner, 1.5T GE 
TwinSpeed scanner and a 3.0T GE VH/i scanner. A standard extremity coil was used for 
imaging the knee and ankle at 1.5T, a prototype phased array multi-coil was used for the 
knee and ankle at 3.0T, and a standard phased array torso coil were used for imaging the 
heart, pelvis and abdomen at 1.5T. The study was approved by our institutional review 
board. Prior to imaging, informed consent was obtained from all volunteers. Product 
automated shim routines were used for all imaging. 

3D-SSFP imaging was performed in the knee, ankle and pelvis. Typical imaging 
parameter for the knee and ankle included: N x =256 (fractional readout), N y =192, NSA=1, 
FOV=16cm, slice=1.5mm, N z =32-64, TR=5.6-6.2, and three or four echoes spaced by 
approximately 1ms (1.5T) and 0.5ms (3.0T). Bandwidth at 1.5T was ±125kHz and ±42kHz 
at 3.0T. Imaging parameters in the pelvis included: N x =512, N y =X, NSA=1, FOV=32cm, 
slice=4mm, N z =32, TR=5.5ms, and three echoes (TE=0.9,1. 9,2.9ms). 

To demonstrate the ability of the "Dixon" method to separate fat and water in the 
presence of a heteroegeous magnetic field, 3D spoiled gradient echo (SPGR) images with 
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fat saturation were obtained in the knees and ankles at both 1.5T and 3.0T. Matrix size, 
field of view and slice thickness were the same as SSFP images. Other parameters 
included: TR=50ms, TE=5ms (full echo), tip angle = 40°, and BW = ±16kHz. These 
parameters are based on established reports using fat saturated SPGR imaging of articular 
cartilage [18, 19]. Total scan time for a 64 slice FS-SPGR set of images was 12:04 min. 
Comparison fat-saturated FSE images were also obtained in the abdomen, with the same 
imaging parameters used for the "Dixon" acquisition, except that 3 signal averages were 
used to obtain comparable SNR. 

Cardiac images were acquired at 1.5T with a modified 2D retrospectively gated 
SSFP sequence with the following imaging parameters: N x =224 (fractional readout), 
N y =128, NSA=1, FOV=32cm, slice=8mm, TR=5.2, and three echoes (TE=0.9, 1.9,2.9). 
Twenty phases were acquired through systole and diasotole with a segmentation factor of 
16 for a time resolution of 83ms. Total breath-hold time for one slice was approximately 
20-24s, depending on heart-rate. 

FSE images were acquired in the knee and abdomen. In the knee, imaging 
parameters included: BW=±16kHz, N x =320, N y =224, NSA=1, FOV=16cm, slice=2mm, 
TR=5000ms, effective TE=10ms, echo-spacing=13.6ms and three echo time increments 
(ATE=- 1.0ms, 0ms, 1.0ms). In the abdomen, imaging parameters included: torso phased 
array coil, respiratory triggering, BW=±31.3kHz, N x =384, N y =192, NSA=1, FOV x =34cm, 
FOV y =25.5cm slice=8mm, effective TE=90ms, echo-spacing= 13.6ms and three echo time 
increments (ATE=-lms, 0ms, 1ms). For comparison, fat saturated FSE images were 
obtained with identical imaging parameters and NSA=3. 

SPGR imaging of a silicone-fat-water phantom was peformed at 1.5T using a 
standard head coil. Sylgard 527 Dielectric Silicone Gel (Dow Corning, Midland, MI) was 
obtained and parts A and B were mixed in a ratio of 6:7 and allowed to set in a 15cc and 
50cc vial. Olive oil was placed in a 15cc and 50cc vial and all vials placed in a one liter 
vessel containing tap water. Imaging parameters included: BW=±3 1.3kHz, N x =256, 
N y =256, NSA=4, FOV=16cm, slice=4mm, TR=34ms, TE=4.9, 7.3, 9.7, 12.1ms. 

An off-line reconstruction program written in Matlab 6.0 (Mathworks, Natick, MA) 
was used to perform fast Fourier transform reconstruction of all images. Following 
reconstruction of complex (magnitude and phase) images, estimation of water images and 
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fat images based on the iterative least-squares algorithm described above, was performed. 
Smoothing of field maps was performed with a 3x3 box-car filter. 

Results: 

Examples of source SSFP and FSE images, as well as calculated water and 
calculated fat images are shown in figs. 2-7,8-10. Uniform fat-water separation was 
achieved in all images with relatively short acquisition times. In SSFP imaging, fluid 
appears bright, which is particularly helpful in musculoskeletal imaging, providing an 
arthrographic effect that increases the conspicuity of cartilage defects. For example, a small 
focal cartilage defect in the anterior tibial cartilage of the ankle is easily seen in the water 
image (Fig. 5). 

Cardiac SSFP images demonstrate uniform fat-water separation with good spatial 
and temporal resolution, acquired within one breath-hold (figs. 7). Flow artifact was 
minimal, and excellent fat-water separation was achieved. 

Uniform fat separation was seen with multi-coil SSFP and FSE imaging in pelvis 
and abdomen (figs. 6,10). A rim of high signal is seen anteriorly within the FSE calculated 
water image (fig. 10b). This was caused by a small misregistration between sequential 
source images that were acquired during free breathing with respiratory triggering. 

Excellent separation of silicone-fat-water was also achieved in a phantom using 
SPGR imaging (fig. 1 1). Four source images with echo time increments=2.4ms were used 
to optimize the SNR of this decomposition, as described in eqs. lOa-c and visualized 
explicitly in fig. lb. 

Discussion: 

In this work we have demonstrated the feasibility of using a multi-point "Dixon" 
chemical shift separation technique that uses multiple images acquired at arbitrary echo 
time increments to separate each species with a unique chemical shift. An iterative least- 
squares fitting algorithm was described and implemented with pulse sequences such as 
SSFP and FSE that benefit from short echo time increments to maintain good image 
quality. This allows for optimization of SNR performance as well as tradeoffs between 
acquisition parameters and SNR. In addition, phase-unwrapping algorithms are not 
necessary with this approach; such algorithms are commonly used with "Dixon" fat-water 
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separation techniques and add tremendously to the complexity of the reconstruction 
algorithm. The approach was also extended to multi-coil acquisitions with minimal 
increase in complexity. In vivo examples were shown in the knee, ankle, pelvis, and heart 
demonstrating excellent fat-water separation. In all cases, uniform separation of all 
chemical species was achieved with excellent image quality and high SNR. An added 
benefit of "Dixon" techniques over fat-suppression techniques is that the source images, as 
well as fat images are available and can contribute to the diagnostic value of the study. 

Artifacts will result if source images are not registered with one another, as was 
seen in fig. 10. Registration of source images in musculoskeletal imaging was not 
problematic, and registration in cardiac imaging was achieved by acquiring all images 
within a breath-hold. In the presence of arrhythmias, misregistration of cardiac source 
images could occur. Imaging within the pelvis was also not problematic, although imaging 
within the adbomen was more difficult because of misregistration from respiratory motion. 
Possible solutions include reduced scan times to acquire images within a breath-hold or 
interleaving echoes for different source images to register data more closely. 

The mathematics of this technique was generalized to include systems with multiple 
chemical species such as silicone/fat/water. Simultaneous separation of these three species 
may be helpful in the assessment of breast implant integrity, and an example of 
silicone/fat/water separation was demonstrated in a phantom. Prior work with traditional 
three-point "Dixon" imaging has exploited the fortuitous relative chemical shifts of 
silicone, water, and fat, permitting the separation of water and fat into one image and 
silicone into a second [17]. 

The noise analysis presented above assumed that the field map is known precisely. 
In general, this is not the case, and the estimate of the field map itself is calculated from 
source images, which contain noise. This in turn will degrade the noise performance of the 
final estimates of water and fat images, and a complete description of the noise 
performance of multi-point fat- water decomposition must reflect the fact that the field map 
is an unknown quantity. The effect of estimating the field map degrades the noise 
performance of the water and fat estimation from that described in eq. 9. 

Other factors in a noise analysis must also be considered. For example, if it 
assumed that the field map is smoothly varying in space, then spatial smoothing of the final 
estimate of the field map will improve the SNR performance of the estimation algorithm. 
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An alternative to smoothing the calculated field map is to smooth the source images, which 
is an approach that has been described previously [20, 21]. In addition, the combination of 
separate field maps calculated in a multi-coil algorithm further improves SNR performance 
of field map estimation by combining multiple contributions to the local field map. The 
description above, summarized in the special case described in eq. 9 is very useful, 
however, as it acts as an upper limit for the best achievable SNR performance for the 
estimation of water and fat, as well as describing the effects of echo time increment. 
Improvements in the noise performance of the field map estimation through smoothing and 
multi-coil acquisitions will improve the NSA for water and fat, but cannot exceed this 
upper limit. A detailed mathematical description that includes the effects of field map 
estimation, smoothing of field maps or source images, and multi-coil acquisitions is 
complex and beyond the scope of this paper. 

The extension of the least-squares iterative algorithm to multi-coil acquisitions was 
straightforward. The calculated field map for each coil is independent of local phase shifts 
that may differ between coils, allowing easy combination of all field maps through well- 
described multi-coil combination methods. This obviated the need for complex phase 
unwrapping algorithms applied to the multi-coil source data. Combination of field maps 
using multi-coil acquisitions also improves the local SNR of the field map, improving noise 
performance of water and fat estimation. 

The use of shorter echo time increments described in this work helps to alleviate the 
large increases in TR with SSFP and echo spacing with FSE, that are necessary with 
standard [5] and modifed [12] three-point "Dixon" techniques. In general, any technique 
that shifts echo times will increase TR and echo spacing and does limit this approach in 
circumstances where the TR (SSFP) and echo spacing (FSE) must be kept very short. 
Recent work by Ma et al describes a three point "Dixon" method that acquires FSE source 
images at different effective echo times through the use of shifted fractional echoes and a 
reconstruction algorithm akin to homodyne reconstruction [22, 23]. This allows acquisition 
of data with echoes that are sufficiently spaced to separate fat and water, with minimal 
increases in the echo spacing. The same approach should be effective with the work 
described above, although the increase in echo spacing was relatively minimal with the low 
bandwidths that were used, and this approach would be most beneficial in FSE applications 
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that use high bandwidths. Application to SSFP imaging would help maintain short TR's, 
preventing image degradation from banding caused by field heterogeneities. 

Although fat-water separation with the method described in this work is efficient 
from an SNR perspective, acquisition of a minimum of three echoes can be limiting, 
especially in cardiac imaging where all data must be acquired within a breath-hold. There 
are many strategies that can be employed to reduce total acquisition time. Many of these 
approaches assume that the field heterogeneity map varies smoothly varying across the 
object, and sampling of the central portions of k-space may be adequate to estimate the 
field map adequately. A simple method for scan time reduction would be the acquisition of 
some source images with a reduced matrix size. For example, in a three-echo acquisition 
scheme, reduction in the number of k y (phase-encoding) lines of one source image would 
permit calculation of a low resolution field heterogeneity map, while reducing total scan 
time. With this low-resolution field map and two full-resolution source images, estimates 
of water and fat images could then be made. 

Additional time savings could be made with pulse sequences that acquire multiple 
images in a time sequence, such as cardiac CINE imaging, in which the field maps can be 
assumed to vary little in time [24]. Field maps that are calculated for one image could be 
used for fat-water separation in multiple other images, reducing the amount of required 
data dramatically. 

The acquisition of multiple echoes within one TR, such as a multi-echo SSFP 
sequence [25] or gradient-echo spin-echo (GRASE) [26] could also exploit the evolving 
phase of spins with different chemical shifts. The mathematics of such reconstruction 
schemes are more complex and must consider additional phase shifts and differential time 
delays between echoes acquired with different gradient polarities [27]. 

In systems where low bandwidths are used and significant chemical shift artifact is 
present in the readout direction, "Dixon" imaging could potentially be used to correct for 
chemical shift artifact. This could be performed by decomposing the chemical species as 
described in this work, then each calculated image would be shifted an amount depending 
on the it's chemical shift and the acquisition bandwidth. Correction for distortion from 
susceptibility could also be performed because the field map is known. Finally the 
corrected images would be recombed into a corrected image free from chemical shift and 
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distortion in the readout direction. A similar approach using subsequent acquisition of 
selective water and fat images using spectral-spatial pulses has been described [28]. 

Conclusions: 

The iterative least-squares multi-point fat-water separation method presented in this 
work describes a new approach that permits the use of variable echo time increments that 
can be exploited to optimize SNR or improve speed performance for demanding sequences 
such as SSFP and FSE. The approach can be used to separate systems with three distinct 
chemical shifts such as silicone, fat, and water, and can also be used with multi-coil 
imaging. 
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Appendix A 



For n=l, .... N, eq. 4 can be written in matrix format, 
S = Ap 

where, 

S = [s 



A = 





... S N Oj 










Pi 


Pi Pi 




AY. 










C 2\ 


-d 2l 








-d\ 2 


C 22 


-d 22 


C M2 


~ d M2 


C \N 


-d lN 


C 2N 


-d 2N 


C MN 


~ d MN 


dn 


c n 






d M\ 


C M\ 


d n 




d 2 2 




d M2 


C M2 


d\N 




d 2N 


C 2N 


d MN 


C MN _ 



(A.1) 



(A.2) 



where rows 1 to N of matrix A are used to calculate the real components of the signal, and 
rows N+l to 2JV are used to calculate the imaginary components of the signal. 

Using a well-described least-squares fitting approach for linear systems of 
equations [29], it can be shown, 



p = (A r A)-A'S 



(A.3) 



where p = [p* p[ Pi f>i- P* PmY » determining initial estimates of each 
chemical species. 
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Appendix B 



If = p* + Ap* and /?; = f>j + Apj (/=i M), and ^ = y/ Q + Ay/ , then eq. 2 can be 

written, 



f M 



(B.l) 



Dividing each side by e /27r ^'« , and using the Taylor approximation e i2 * Av " n « 1 + /2;rA^, J , 
such that, 

Rearranging eq. B.2, and splitting into real, 
and imaginary components, 

+ #0= 2 ^t (# c * -p'jd Jn ) + ±(^d JK + A P ' jCjn ) (b.4) 

where i* and !„' are defined in eqs. B.3 and B.4. Arranging in matrix format for n-\,..., N 

S « By ( B - 5 ) 
where! = [?,* If ... s R N $[ H - > 

y = [A^ Apf Ap[ bp* Api ... Ap R M Ap'J , g% = 2* r„ -p'jC Jn ) 



and g' jn = 2nt n Y i {p*c J „- P ' J d J X such that, 



B = 



ft - , 
ft* 2 



a* 

ft', 
ft' 2 



•-11 
C,2 



•-22 
C 2N 

d 2 \ 
d 22 



-d 2 
-d 



22 



-d 



IN 
C 21 
C 2 2 



'A/2 



'A/1 



A/2 



-Ml 



-Ml 



Sw d \ 



IN 



d MN C MN J 



(B.6) 
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For M =y jv, e q. B.5 is also a linear system of equations, and estimates of y can be 

calculated [29], 



y = (B r B) _l B r S 
which is used to determine Ay/, Ap R . , and Ap'j 



(B.7) 



Appendix C 

In the special case of only water and fat (A^= relative fat-water chemical shift) 
and the receive/transmit frequency of the scanner is set to the water resonance, matrices A 
and B become, 



and 



A = 



1 
1 

1 
0 
0 



B 



g. 

gl 



g« N 



g\ 



g's 



0 
0 

0 

1 
1 



1 
1 

1 

0 
0 



df 
d? 



0 
0 

0 

1 
1 



-d? 



d S Z 



-df 
-d? 



cr 



C.l 



where cf=co S (2^J, d? = sin(2^J, g* = 2*/.(- -fy.) and 
g' n =2xt„{p* w +p* f c n -p l f d n ) are the matrix elements. 
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TE Increment (ms) 



Figure 1: a) Dependence of effective number of signal averages (NSA) on echo time 
increment (ms) from fat-water chemical shift for three-echo (solid) and four-echo (dashed) 
acquisitions at 1.5T. b) NSA for water, fat, and silicone for a four-echo acquisition to 
resolve all three species at 1.5T. 
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Figure 2- Sagittal SSFP images obtained in the knee of a human volunteer at 1.5T. Source 
(a), calculated water (b), calculated fat (c) images are shown in comparison to a fat- 
saturated SPGR image (d). Excellent fat-water suppression was obtained with the bbtv 
images, and joint fluid appears bright (arrows), providing an arthrography effect 
potentially improving the conspicuity of cartilage defects. Four source images (TE- 1.2 
2.1, 3.0, 3.9ms) were acquired and TR=6.1ms. Image matrix size was 256x192x64, and 
bandwidth was + 125kHz. 
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Figure 3- Sagittal SSFP images obtained in the ankle of a human volunteer at 3.0T. Source 
(a) calculated water (b), calculated fat (c) images are shown in comparison to a fat- 
saturated SPGR image (d). Compared to FS-SPGR images, excellent fat-water suppression 
was obtained with the SSFP images, and joint fluid appears bright. Three , source images 
were acquired (TE=1.4, 1.9, 2.4ms) and TR=5.6. Image matrix S1 ze was 256x192x32, and 
bandwidth was ±42kHz. 
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Figure 4- Sagittal SSFP images of the knee of a human volunteer at 3.0T with a) source b) 
calculated water, and c) calculated fat images. Uniform fat-water separation * identified 
and joint fluid appears bright. Three source images were acquired (TE-1.4, 1.9, 2.4ms;, 
and TR=5.6ms. Image matrix size was 256x192x32, and bandwidth was ±42kHz. 
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Figure 5: Sagittal SSFP images of an ankle at 1.5T with a) source, b) calculated water, and 
c) calculated fat images. Uniform fat-water separation is identified and joint fluid appears 
bright. A small defect in anterior tibial cartilage (arrow) is seen and well delineated by joint 
fluid which appears bright. Four source images were acquired (TE=1.2, 2.1, 3.0, 3.9ms), 
and TR=6.1ms. Image matrix size was 256x192x64, and bandwidth was ±125kHz. 
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a a -vtoi <j<;fp imaees obtained in the pelvis of a human volunteer at 1.5T. Source 

£££ and good Lige Tl». ^T^T^ttSSS i£ 

2.8ms) and TR=5.4ms. Image matrix size was 512 x 250 x sl, ana o 



±125kHz. 
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Figure 7: Short axis CINE SSFP cardiac images acquired with torso phased array coil. 
Source image (left), calculated water image (middle) and calculated fat image (right) are 
shown at one (of 20) end diastolic phase. Uniform fat-water separation was consistently 
achieved in all slices and phases. Three source images were acquired (TE=0.9 , 1 .9, 2.9ms) 
and TR=5.2ms. Image matrix size was 224x128, and bandwidth was ±125kHz. 
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Figure 8: Sagittal proton density FSE images of a knee of a human volunteer acquired at 
1.5T with a) source, b) calculated water, and c) calculated fat images. Uniform fat-water 
separation is achieved and joint fluid appears bright. Three source images were acquired 
(shift in TE = -1.0, 0, 1.0ms), and TR=5000ms. Effective TE was 10ms and image matrix 
size is 320x224. Echo spacing was 13.7ms, and bandwidth was ±20. 8kHz. 
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Figure 9: Axial T 2 -weighted density FSE images of a knee of a human volunteer acquired 
at 1 5T with a) source, b) calculated water, and c) calculated fat images. Uniform fat-water 
separation is achieved and joint fluid appears bright. Three source images were acquired 
(shift in TE = -1.0, 0, 1.0ms), and TR=5000ms. Effective TE was 70ms and image matrix 
size is 384x192. Echo spacing was 17.5ms, and bandwidth was ±15.6kHz. 
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Figure 10: Axial T 2 -weighted fast spin-echo images obtained from the abdomen of a human 
volunteer at the level of the spleen at 1.5 with a) source, b) calculated water, c) calculated 
fat images. Images were acquired with the torso phased array coil during free breathing 
using respiratory triggering. Uniform fat-water separation was achieved everywhere except 
at the anterior skin surface where misregistration between subsequent scans caused 
calculation errors (arrows). Three source images were acquired (shift in TE=- 1.0,0, 1.0ms). 
Effective TE was 90ms, and image matrix is 384x192 with a 3 / 4 FOV. Echo spacing was 
1 1.9ms and bandwidth was ±3 1.3kHz. Fat-saturated T 2 -weighted fast spin-echo image with 
the same imaging parameters and NSA=3 is also shown (d). 
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Figure 1 1 : SPGR a) source, b) calculated water, c) calculated fat, d) calculated silicone 
images in a phantom. Four images were acquired (TE=4.9, 7.3, 9.7, 12.1ms) and TR=34ms. 
Other parameters included: 10° flip angle, 256x256 matrix, 16cm field of view, ±31. 3kHz 
bandwidth, and 4mm slice thickness. 
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